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The ground state of the stacked triangular antiferromagnet CsFeBr3 is a spin singlet due to 
the large single ion anisotropy D{S Z ) 2 . The field-induced magnetic ordering in this compound 
was investigated by the magnetic susceptibility, the magnetization process and specific heat 
measurements for an external field parallel to the c-axis. Unexpectedly, two phase transitions 
were observed in the magnetic field H higher than 3 T. The phase diagram for temperature 
versus magnetic field was obtained. The mechanism leading to the successive phase transitions 
is discussed. 
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§1. Introduction 

The spin frustration effect often plays an important 
role in the magnetic ordering process and magnetic ex- 
citations. In the hexagonal antiferromagnets of ABX3- 
type with the CsNiCi3 structure, magnetic B 2+ ions form 
infinite chains along the c-axis and triangular lattices in 
the basal c-plane. Since the exchange interaction in the 
c-plane is antiferromagnetic, they behave as triangular 
antiferromagnets (TAF) at low temperatures. Because 
of the spin frustration effect being characteristic of TAF, 
together with the quantum effect, a rich variety of phase 
transitions have been observed in the hexagonal ABX3 
antiferromagnets. ^ 

Low-temperature magnetic properties of AFeX3 sys- 
tems are described by the pseudo spin S — 1 anisotropic 
XX Z model with the large easy-plane anisotropy of the 
form D{S Z ) 2 due to the crystalline field. 2 ) In RbFeCl 3 , 
exchange interactions overcome the anisotropy, so that 
RbFeCl3 undergoes magnetic phase transition in the ab- 
sence of magnetic field. On the other hand, CsFeCl3 has 
a singlet ground state at zero field, because exchange in- 
teractions are not sufficiently strong to produce the mag- 
netic ordering. The ordering process in RbFeCl3 is not 
simple. The phase transition occurs from the paramag- 
netic state to the commensurate (C) ground state with 
the 120° structure through two different incommensurate 
(IC) states. 3 ' 4 ' The IC-C phase transition occurs due to 
the competition between the antiferromagnetic exchange 
interaction in the c-plane and the dipolc-dipolc (D-D) 
interaction, the latter of which is enhanced by the fer- 
romagnetic exchange interaction along the c-axis. 5 ' 6 - ) A 
similar IC-C phase transition was observed in CsFeCl3 



Table I. Interaction parameters for CsFeBr3 in the unit of K. 
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when magnetic fields H are applied parallel to the c- 
axis. 3 ' 7 ' 8 -' In addition, in RbFeCi3, the quantum fluctu- 
ation produces successive phase transitions in magnetic 
fields perpendicular to the c-axis. 9_11 ^CsFeBr3 is isostruc- 
tural with RbFeCl 3 and CsFeCl 3 . 12) However, CsFeBr 3 
differs from them in terms of the exchange interaction 
along the c-axis. CsFeBr3 has an antiferromagnetic in- 
trachain interaction. Thus, the D-D interaction should 
be negligible. As observed in CsFeChj, the ground state 
of CsFeBr 3 is a spin singlet due to the large single ion 
anisotropy. 13 ' 14 -* The effective Hamiltonian of CsFcBi-3 
at low temperatures for H \\ c may be written as 

H = J2[D(sn 2 -gmHS*] 

i 

chain 
plane 

+ E [iHstsr + s-s+) + 2jfs;s;] , 

(1.1) 

where Sf (a — x, y, z) is the spin-1 operator on the i-th 
Fe 2+ site, and jj^ and jf' X are the antiferromagnetic 
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Fig. 1. The temperature dependence of the magnetic susceptibil- 
ity M/H in CsFcBr 3 at 1 and 7 T for H || c. 



exchange interactions along the c-axis and in the c-planc, 
respectively. 

Magnetic excitations in CsFeBr3 have been extensively 
investigated by means of neutron inelastic scattering ex- 
periments at zero and finite magnetic fields. 14 ~ 19 ) The 
lowest excitation occurs at Q — (1/3,1/3,1) and its 
equivalent reciprocal points with the excitation energy 
of A = 0.11 THz. 18 ) The results were analyzed on the 
basis of the Hamiltonian of eq. (1.1), and the interaction 
parameters were determined, as listed in Table I. 

When an external field H is applied along the c- 
axis, CsFcBr 3 undergoes a magnetic phase transition for 
H > 4 T. 16 ) The magnetic structure in the ordered phase 
is a triangular structure characterized by an ordering 
vector Q — (1/3, 1/3, 1), which was at first assumed to 
be the 120° spin structure in the basal plane. 16 ' ^ The 
stacked TAF with XY spins exhibits a phase transition 
from the paramagnetic state to the ordered state with the 
120° spin structure. Recently, the nature of this phase 
transition has been under controversy. Kawamura 21 ' 22 ' 
argued that the phase transition is of the second order, 
and belongs to a new chiral XY universality class. His 
theory was supported by neutron scattering and specific 
heat experiments performed on CsMnBr3. 23 ~ 26 - ) How- 
ever, it was shown later that the phase transition must 
be of the first order. 27 ' 28 ) This problem motivated us to 
study the field-induced magnetic ordering in CsFeBr 3 . 
Since CsFeBr 3 has a large easy-plane anisotropy, which 
overcomes all exchange interactions, the field-induced 
phase transition for H \\ c may have the same nature 
as that of the XY stacked TAF. In order to investigate 
the field-induced magnetic phase transition in CsFeBr3, 
we have performed magnetization and specific heat mea- 
surements. As shown below, unexpected two-step phase 
transitions were observed. 

The arrangement of this paper is as follows: In §2, 
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Fig. 2. The low-temperature magnetic susceptibility in CsFeBr3 
at various external fields for H II c. 



the experimental procedures are described. The exper- 
imental results and their analyses are presented in §3. 
Sections 4 and 5 are devoted to the discussion and con- 
clusion, respectively. 

§2. Experimental Procedures 

Single crystals of CsFeBr3 were grown by the vertical 
Bridgman method from the melt of an equimolar mix- 
ture of CsBr and FeBr2 sealed in evacuated quartz tubes. 
After weighing, they were placed in a quartz tube and 
dehydrated by heating in vacuum at ~150°C for three 
days. The temperature at the center of the furnace was 
set at 650°C, and the lowering rate was 3 mm-h -1 . Sin- 
gle crystals of 1^5 cm 3 were obtained. The crystals ob- 
tained were identified to be CsFeBr3 by X-ray powder 
diffraction. 

The specific heat measurements for single crystal of 
CsFeBr3 were carried out at RIKEN down to 0.6 K 
in magnetic fields up to 10 T using a Mag Lab HC mi- 
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Fig. 4. The phase boundary in CsFeBr3 for H \\ c determined 
from the results of the temperature variations and field varia- 
tions of the magnetization and the specific heat. Solid circles 
and rectangles denote the transition points determined from the 
temperature dependence of susceptibilities and the magnetiza- 
tion curves, respectively. Diamonds and triangles mean the tran- 
sition points determined from the temperature and field varia- 
tions of specific heat, respectively. Solid and dashed lines are 
the fits with eq. (3.1) with </>i = 2.2(1) and H N1 (0) = 2.89(2) 
T, and <j> 2 = 1.8(1) and H N2 (0) = 2.88(2) T, respectively. 



Fig. 3. Magnetization curves in CsFeBr3 at various temperatures 
for H || c. The inset shows dM/dH versus H. 



crocalorimeter (Oxford Instruments) in which the relax- 
ation method was employed. The magnetizations were 
measured down to 1.8 K in magnetic fields up to 7 T 
using a SQUID magnetometer (Quantum Design MPMS 
XL). The high-field magnetization measurement was per- 
formed using an induction method with a multilayer 
pulse magnet at the Ultra-High Magnetic Field Labo- 
ratory, Institute for Solid State Physics, The University 
of Tokyo. Magnetization data were collected at 1.6 K in 
magnetic fields up to 40 T. In these experiments, mag- 
netic fields were applied along the c-axis. 

§3. Experimental Results 

3.1 Magnetic susceptibility and low- field magnetization 
We first measured the magnetic susceptibilities in 
CsFeBr 3 at H = 1 and 7 T for H || c. The results 
are shown in Fig. 1. With decreasing temperature, the 
susceptibility \ (= M/H) exhibits a broad maxima at 
about 15 K and then decreases. The susceptibility for 
H = 1 T decreases monotonically toward zero, while the 
susceptibility for H = 7 T displays a cusplike minimum 



at T — 2.62 K, which indicates the magnetic phase tran- 
sition. Figure 2 shows the low-temperature susceptibility 
in CsFeBr 3 measured at various external fields for H \\ c. 
We assign the temperature at which there is an inflection 
point in the susceptibility as the transition temperature 
Tn{H). The arrows in Fig. 2 indicate the transition tem- 
peratures, which arc slightly lower than the temperatures 
of the susceptibility minima. With increasing external 
field, the transition temperature increases. 

Figure 3 shows the magnetization curves in CsFeBr3 
measured at various temperatures for H \\ c. dM/dH 
versus H is presented in the inset. A sharp peak 
anomaly, which is indicative of the phase transition, can 
be observed in dM/dH versus H . We assign the field at 
which dM/dH exhibits the peak anomaly as the tran- 
sition field Hn{T). The arrows in Fig. 3 indicate the 
transition fields. With increasing temperature, the tran- 
sition field increases. In spite of the singlet ground state, 
the slope of the magnetization curve for H < H^ is fairly 
large even at T = 1.8 K. This is because of the large Van 
Vleck paramagnetism due to the crystalline field. 

In the present magnetization measurements, a single 
phase transition was clearly observed. In Fig. 4, the 
phase transition temperatures T^(H) and fields Hn(T) 
were plotted as solid circles and rectangles, respectively. 
Because the transition points determined from T^(H) 
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Fig. 5. The temperature dependence of the specific heat in 
CsFeBr3 at various magnetic fields for H \\ c. The values of the 
specific heat are shifted upward by 1 J/mol K with increasing 
external field. 



and Hn{T) lie almost on the same line, they are consis- 
tent with each other. 



3.2 Specific heat 

Figure 5 shows the total specific heat C in CsFeBr3 
measured at various external fields for H \\ c. For H = 3 
T, no anomaly is observed down to T — 0.6 K, while for 
H = 3.5 T, the anomaly composed of two peaks can be 
observed. With increasing external field, the two peaks 
shift toward the high-temperature side, and separate in- 
creasingly. They are well defined for H > 7 T. We con- 
firmed that the sample had not been cracked into two 
pieces during the cooling, and that the values of the spe- 
cific heat are reproducible using another sample. Thus, 
the two-peak anomaly in C is not due to a sample prob- 
lem, but is intrinsic to the present system. Based on the 
result of the specific heat measurements, we can con- 
clude that the magnetic ordering in CsFeBr3 occurs in 
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Fig. 6. The field dependence of the specific heat in CsFeBr3 at 
0.62 K for H II c. 



two steps. The arrows in Fig. 5 denote the phase tran- 
sition temperatures Tni and Tn2- 

In order to determine the phase transition points for 
H < 3.5 T, which are difficult to determine from the 
temperature scan, we measured the field dependence of 
C for H | c. Figure 6 shows the measurements at 
T = 0.62 K. A shoulder and cusplike anomalies are vis- 
ible at H m = 3.05 T and H m = 3.25 T.In Fig. 4, the 
phase transition points determined from the specific heat 
measurements are plotted as diamonds and triangles. It 
is evident that the outer phase boundary determined 
from Tni and H^i coincides with the phase boundary 
determined from the magnetization measurements. The 
temperature interval of the intermediate phase increases 
with the magnetic field. We label the intermediate phase 
and the low-temperature phase inside the inner bound- 
ary phases I and II, respectively. When the magnetic 
field is slightly higher than the gap field H g correspond- 
ing to the excitation gap, both phase boundaries may be 
represented by the power law 



[H Na (T) - H Na (0)] ocT*- 



(3.1) 



with a =1 and 2. We fit eq. (3.1) to the data for T < 2 
K, which is lower than half the gap temperature A/fcs ~ 
5.3 K. 18 ' Solid and dashed lines in Fig. 4 are fits with 
fa = 2.2(1) and 7T N i(0) = 2.89(2) T, and fa = 1.8(1) 
and i?N2(0) = 2.88(2) T, respectively. Since ff N i(0) « 
^N2(0), we can conclude that the two phase boundaries 
meet at zero temperature. Thus, the gap field H g for 
T = is estimated to be H g = H m (0) = 2.89 T. 

3.3 High- field magnetization 

In order to determine the phase boundaries in high 
fields, we performed high-field magnetization measure- 
ments. Figure 7 shows the magnetization curve and 
dM/dH versus H measured at T = 1.6 K. Sharp anoma- 
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lies indicative of phase transitions are observed at H c = 
3.54 T and H s = 27.8 T. The transition field H c should 
correspond to H^i observed through low-field magneti- 
zation and specific heat measurements. H s is the satu- 
ration field. 

The magnetization curve for H > H s has a finite 
slope due to the large Van Vleck paramagnetism, i.e., 
Xvv = 4.64xl0~ 3 emu/mol. The saturation magneti- 
zation is evaluated to be M s = 2.65 ± 0.10 /^B/Fe 2+ 
by extrapolating the magnetization curve for H > H s 
to zero field. The value of M s is almost the same as 
M s = 2.6 /ib /Fe 2+ in CsFeCl 3 . 8 ) From the value of M s , 
the g-factor for H || c is obtained to be g\\ = 2.65 ± 0.10. 
The g-factor can also be evaluated to be g\\ — 2.7 us- 
ing the relation of A = g^H s with A = 0.11 THz 18 ) 
and H g — 2.89 T. Both ^-factors coincide and agree with 
previously obtained g\\ = 2.6. 17 ' 

The transition points obtained from the high-field 
magnetization measurements are plotted as solid trian- 
gles in Fig. 8 together with those obtained from the 
low-field magnetization and specific heat measurements. 
The solid lines are guides for the eyes. The dashed lines 
denote the expected phase boundaries. 

§4. Discussion 

The field-induced magnetic ordering in the singlet 
ground state magnet with the large D term was first 
investigated by Tsuneto and Murao 29 -' on the basis 
of the mean-field theory. They demonstrated that a 
single-ordered phase exists in a closed area in the mag- 
netic field vs temperature (H — T) diagram. Such 
field-induced phase transition was observed in sev- 
eral nickel compounds, e.g., Ni(C 5 H 5 NO) 6 (ClO 4 )2, 30 ' 31) 
Ni(N0 3 ) 2 • 4H 2 32 > and Ni(N0 3 ) 2 ■ 6H 2 0, 33 ) all of 



which are unfrustrated systems. 

Since the intrachain exchange interaction Jo is anti- 
ferromagnetic in CsFeBr3, the dipole-dipole interaction 
should be negligible. Thus, it has been considered that 
CsFeBr3 would undergo a single phase transition from 
the paramagnetic phase to the ordered phase with the 
120° spin structure in the basal plane in high magnetic 
fields. However, the present study revealed that the field- 
induced phase transition occurs in two steps. This is an 
unexpected finding. The phase diagram for H \\ c was 
determined as shown in Figs. 4 and 8. 

The magnetic ordering in CsFeBr3 was investigated by 
means of neutron elastic scattering in magnetic fields up 
to 6 T. 16 ' 34 ) The magnetic Bragg reflections were ob- 
served at Q = (1/3,1/3,1) and (2/3,2/3,1). Thus, it 
is certain that a triangular spin structure is realized in 
the basal plane in the phase II. However, it was also sug- 
gested from the magnetic excitations in phase II that 
the triangular spin structure is not the regular 120° 
structure but is a slightly modified one. 34 - 1 The transi- 
tion fields determined from the field dependence of the 
Bragg reflections for T < 1.6 K appear to be on the 
outer phase boundary in Fig. 4 (or Fig. 8). The dis- 
tinct anomaly that is indicative of the second transi- 
tion was not observed in the field variation of the mag- 
netic Bragg reflections. Since the field interval of the 
intermediate phase I is less than 0.3 T for T < 1.6 K, 
it seems difficult to distinguish phase I. In the closely 
related compound RbFeBr3, the total of the exchange 
interactions overcomes the single ion anisotropy energy 
(see eq. (A. 13) in the Appendix), so that the magnetic 
ordering occurs even at zero field. RbFeBr3 undergoes 
two magnetic phase transitions at Tni = 5.61 K and 
Tn2 = 2.00 K at zero field. 35 ) This compound has two 
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Fig. 8. The phase diagram in CsFeBr3 for H \\ c determined by 
the present measurements. The solid lines are guides for the 
eyes. The dashed lines denote the expected phase boundaries. 



structural phase transitions at T\ = 109 K and Tq = 34.4 
K, which are accompanied with the shift of -FeBr3- 
chains along the c-axis. 36 ' 37 ^ The chemical unit cell be- 
low T\ is described by enlarging that of CSMCI3 to v^o, 
y/3a, c. Consequently, there exist two different inter- 
chain exchange interactions, J\ and J{, which are equiv- 
alent in the undistorted CSN1CI3 structure. The succes- 
sive magnetic phase transitions in RbFeBr3 arise from 
Ji 7^ t/(, 35 ^ which are explainable in the framework of 
the mean-field theory. A collinear and a triangular spin 
structures are realized in the intermediate- and the low- 
temperature phases, respectively. The successive mag- 
netic phase transitions due to the same mechanism were 
also observed in RbVBr3. 38 'Such structural phase transi- 
tion has often been observed in hexagonal RbBBr3 com- 
pounds, e.g., RbFeBr 3 , 36 ' 37 ) RbVBr 3 , 39 ) RbMnBr 3 40 ) 
and RbCoBr3, 41 ' but not in CsBBr3 compounds ex- 
cept for the Jahn- Teller system CsCrBr3. 42 ' CsFeBr3 
has the same crystal structure as CsNiCl3 (space-group 
P63 / mmc) at room temperature, 12 '' which we confirmed 
by X-ray powder diffraction. In order to check whether 
there is a structural phase transition in CsFeBr3, we mea- 
sured the dielectric constant between 4.2 K and 300 K at 
zero magnetic field. However, no anomaly indicative of a 
structural phase transition was detected. Thus, the crys- 
tal structure remain highly symmetric down to helium 
temperatures. Therefore, the successive field-induced 
magnetic phase transitions in CsFeBr3 should not be at- 
tributed to Ji 7^ J[ due to the structural phase transi- 
tion. However, there is a possibility that the magnetoe- 
lastic coupling produces the situation J\ 7^ J[ below Tni 



Fig. 9. Magnetization curve for H \\ c calculated by the mean- 
field approximation with gn = 2.65 and the magnetic parameters 
obtained by Dorner et aZ. 14 ' (see Table I). The solid and dashed 
lines are the results in the cases when the highest — 1 ) state is 
included and excluded, respectively. The thick solid line is the 
magnetization curve observed in CsFeBr3 at T = 1.6 K. 



in magnetic fields. 

The two-step transitions observed in CsFeBr 3 cannot 
be described within the framework of the mean-field the- 
ory for the undistorted crystal structure, because it pre- 
dicts a single phase transition. In the triangular antifer- 
romagnet with the singlet ground state as CsFeBr3, both 
the spin frustration and the quantum fluctuation may 
play important roles in the field-induced magnetic phase 
transition. Thus, not only the transverse spin fluctua- 
tion, but also the longitudinal spin fluctuation, which can 
be neglected in the classical spin system, will be respon- 
sible for the ordering process. We infer that the inter- 
play of the spin frustration and the quantum fluctuation 
produce the two-step phase transitions in CsFeBr3. Mi- 
croscopic measurements as neutron scattering and NMR 
are necessary to elucidate the nature of the phase tran- 
sitions. 

Finally we discuss the magnetization process for H || c. 
Assuming that the spin structure in the ordered state is 
an umbrella structure, in which the ^-components of 
spins form the 120° structure in the basal plane, we can 
calculate the magnetization curve for T — by the mean- 
field approximation shown in the Appendix. The thin 
solid line in Fig. 9 is the magnetization curve calculated 
with <7|| = 2.65 and the magnetic parameters obtained by 
Dorner et al. 14 ' (see Table I). It is noted that the ground 
state at zero field is gapless for the magnetic parameters 
obtained by Visser and Harrison. 15, 19 ) The dashed line is 
the M — H curve obtained by neglecting the highest | — 1 ) 
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state. In Fig. 9, we also plotted the experimental result 
obtained at T = 1.6 K, from which the magnetization 
due to the Van Vleck paramagnetism was subtracted. 

The saturation field H s and the gap field H g are ob- 
tained to be H s = 33.7 T and H g = 1.8 T, which agree 
roughly with the experimental saturation field and the 
gap field. In more detail the calculated saturation field 
H s is somewhat larger than the experimental value of 
H s = 28 T, which is expected for T — > 0, while the 
calculated gap field H g is somewhat smaller than the 
experimental value of H g = 2.89 T. The exchange pa- 
rameters should be reduced slightly against the single 
ion anisotropy D to reproduce the experimental satura- 
tion field and the gap field within the framework of the 
mean-field approximation. 

The calculated M — H curve is the convex function 
for H g < H < H s and exhibits remarkable rounding in 
the low field region. This behavior is ascribed to the 
contribution of the | — 1 ) state. On the other hand, the 
magnetization curve observed in CsFeBr3 is almost linear 
for H c i < H < 21 T, and its slope increases steeply for 
H > 21 T (see also dM/dH versus H in Fig. 7(b)). This 
behavior is characteristic of the S = 1 one-dimensional 
antiferromagnetic system. 43 ' 44 ) The experimental result 
indicates that the magnetization process in CsFeBr3 is 
dominated by the spin correlation in the linear chain 
along the c-axis. 

§5. Conclusions 

We have presented the results of magnetization and 
specific heat measurements on a pseudo spin-1 triangu- 
lar antiferromagnet CsFeBr 3 with a singlet ground state 
due to the large single ion anisotropy. The specific heat 
measurements revealed that the field-induced magnetic 
ordering for H || c occurs in two steps. The phase dia- 
gram was determined, as shown in Figs. 4 and 8. The 
successive phase transitions cannot be understood within 
the framework of the mean-field theory. We suggest that 
the phase transitions are attributed to the interplay of 
the spin frustration and the quantum fluctuation orig- 
inated in the singlet ground state triangular antiferro- 
magnet. 
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Appendix: Magnetization Curve for H \\ c at 

T = 

Using the mean-field approximation, we calculate the 
magnetization curve at T = on the basis of the six- 
sublattice model. We assume that an umbrella spin 
structure is realized in the ordered state, for which all 
the sublattice spins are equivalent. We take the sublat- 
tices 1 ~ 3 in a c-plane and 4 <~ 6 in the next c-plane. 



We describe the basis state for the n-th sublattice spin 

as 

ip n = |O)cos0+(| l)cos0 + | - l)sin0)sin0eTr i(n_1) , 

(A.l) 

for n = 1 <~ 3, and 

i\) n = |O)cos0-(|l)cos0 + | - I)sin0)sin0e 4 r t ("- 1 ) , 

(A.2) 

for n = 4 <~ 6, where | 0) and | ± 1) denote the spin 
states for S z — and ±1, respectively. Angles 9 and 
were introduced to satisfy the normalization condition. 
The average values of spin operators are given by 

{S z n ) = sin 2 9 cos 20, {{S z n f) = sin 2 9 , (A.3) 



{St) = (Sf>* = -(St) = -<5 4 ">* 

= V2 cos 9 sin 9 (cos + sin 0) , (A. 4) 

(si) = (s^r = -(s+) = -(s^r 

= V2 cos 6 sin 9(c~^ 1 cos (j) + e^ 1 sin <j>) , (A. 5) 

(si) - <s 3 -T = -<s 6 + > = -<s 6 ->* 

= V2 cos 6 sin^e^cos + e~^ 1 sin <j>) . (A. 6) 
Therefore, the energy per site E is expressed by 
E = _Dsin 2 8 — g\\^H sin 2 9cos2(j> 

- 2 (2Ji + 3Ji) cos 2 6 sin 2 9(1 + sin 20) 

+ 2 ( 4 + Mi) sin4 c °s 2 20 . (A.7) 

Angles 9 and (p can be determined by dEjdQ = and 
dE/d(j) — 0, which lead to 

sin 2 9 = 

9\\HbH cos 2(f) - D + 2 (2j£ + 3Jf) (1 + sin 2(f)) 



(2Jq" + 3Ji) (1 + sin 2(f) + (4 + 3J|) cos 2 20 



(A.8) 



and 

9\\(J,bH sin 20 - 2 (2 Jq" + 3jj~) cos 2 9 cos 20 

-2 (jj +3 4) sin 2 sin 40 = , (A.9) 

for the ordered state. The magnetization curve is ob- 
tained by solving eqs. (A.8) and (A.9) self-consistcntly. 
In the saturated state, sin 6* = 1 and sin0 = 0. Thus, 
the saturation field H s is given by 

.911 n B H s = D + 2 (2 Ji + 3Ji + 2 4 +6 jf) . (A.10) 

At the gap field H g , sin 6 = 0. Substituting this condi- 
tion into eqs. (A.8) and (A.9), we obtain 



g^ B H g = JD* - AD (2Ji + 3Ji) . (A.ll) 



8 



Yosuke Tanaka et al. 



If we neglect the highest 
we obtain 



1 ) state, i.e., sin0 = 0, the 



<7||^Btfg = D - 2 (2J X + 3Jf) . 



(A.12) 



From eq. (A. 11), the critical value of D, which separates 
the ground state natures in the absence of the magnetic 
field, is given by 



^ = 4(2^ + 37^) . 



(A.13) 



The ground state can have a magnetic ordering for D < 
D c while it is a spin singlet for D > D c . 
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